Low-lying Dirac eigenmodes and monopoles in 4D compact QED 
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Abstract 



We study the properties of low-lying Dirac modes in quenched compact QED at /3=0.99, 
1.01 and 1.03, employing 12 3 x 12 lattices and the overlap formalism for the fermion action. 
We pay special attention to the spatial distributions of the low-lying Dirac modes. Near-zero 
modes are found to have universal anti-correlations with monopole currents below/above the 
critical j3. We also study the nearest-neighbor level spacing distribution of Dirac eigenvalues 
and find a signal of a Wigner-Poisson transition. We make a possible speculation on the 
chiral phase transition in 4D compact QED. 
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£1 ■ 1 Introduction 

CN 

■ The chiral symmetry breaking is one of the most famous nonperturbative phenomena in 
QCD, which have been attracting a great deal of interest for a long time. Chiral symmetry, 

■ which is an approximate global symmetry in QCD, is spontaneously broken by nonperturba- 
tive dynamics of QCD. The broken chiral symmetry leads to large constituent quark masses 
of a few hundred MeV PQ [5J |3], which are responsible for about 99% of mass in the world, 
aside from some unknown factors such as dark matters. 

The spontaneous chiral symmetry breaking is caused by non-zero chiral condensate {tpi/i) > 
which is directly related to non- vanishing spectral density at the spectral origin of the Dirac 
43 ' operator, via the Banks-Cashcr relation 4J. The level dynamics of Dirac eigenvalues is then 

essential for the chiral phase transition. As for QCD, the level statistics of the low-lying 
eigenvalues of the Dirac operator has been studied with much effort and it is known to be 
reproduced by the random matrix theory (RMT) with the same global symmetry. There 
however remain several issues to be studied from the microscopic viewpoint. 

The spontaneous chiral symmetry breaking can be also seen in compact QED in 1+3 
dimensions. Compact QED in a strong-coupling region exhibits the charge confinement as 
well as the spontaneous chiral symmetry breaking like QCD. Its nature has been also studied 
by many groups so far © E M EE EQ E3 EH HI HI EE OH EE EE [IB US HI 123 121 
[251 HE W\- 

The key ingredient for the nonperturbative features in compact QED would be monopoles 
in a microscopic sense. Monopoles play essential roles in the charge confinement and its phase 
transition, and would play important roles also in the chiral symmetry breaking. As a matter 
of fact, neither the confinement nor the chiral symmetry breaking occurs without monopoles' 
degrees of freedom. In this paper, we mainly pay attention to the spatial distributions of 
low- lying Dirac eigenmodes in quenched compact QED, because the low-lying eigenvalue 
dynamics, which is responsible for the chiral symmetry breaking, is considered to have a close 
connection with the low-lying eigcnfunctions. We investigate the spatial correlations between 
low- lying eigenfunctions and monopoles, the key player in the nonperturbative dynamics in 
compact QED. 

The paper is organized as follows. We give the formalism employed in the present analysis 
in Sec. In Sec. El the fundamental features of the low-lying Dirac modes are summed up. 
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We show the lattice QED results in Sec. [U Sec. [S] is devoted to the discussions and the 
speculations based on the present lattice QED results. We finally make a summary in Sec. [SI 

2 Formalism 

2.1 compact QED 

The Wilson gauge action for compact QED is written as 

X fl,V 

with link variables U^(x) = e %6 "^ £ £7(1) and plaquettes P^ lJ (x) = U f _ l (x)U u (x + p,)U^(x + 
v)Ul(x). The constant = corresponds to the coupling constant. Here, the angle O^x) 
ranges from — 7r to n. Defining a plaquette angle 6^ u (x) = 9^{x) + 6 v {x + jl) — 9^{x + 0) — 
y {x) £ (— 47r, 47r], S'qed is represented as 

5qed = /3^^(1-cos^(x)). (2) 

X 

This compact formulation leads to several nonperturbative phenomena, such as the chiral 
symmetry breaking or the charge confinement at small 0. Indeed, at zero temperature, 
compact QED has two phases; the confinement phase and the Coulomb phase. These two 
phases are separated with the critical value C — 1.0111331(21) [24]. The system is in the 
confined phase at < C , while it's in the Coulomb phase at > C . The phase transition 
is weak first order. 

In this paper, we employ 12 3 x 12 lattices at 0=0.99, 1.01 and 1.03. We generate and 
investigate independent 48 gauge configurations at each 0, which are generated with the 
standard Wilson gauge action imposing the periodic boundary conditions in all the directions. 
The system at 0=0.99 (1.03) clearly lies in the confinement (Coulomb) phase, respectively, 
but the 12 3 x 12 system at 0=1.01 is marginal. "Finite temperature" phase transition in 
compact QED was extensively investigated in Ref. [26], and the "transition temperature" 
l/T at 0=1.01 is found to be about l/T ~ 6. Then the 12 3 x 12 system at 0=1.01 is 
considered to be in the confinement phase. 

As we mentioned in Sec. [T] the key ingredient in compact QED in the confinement phase 
is monopoles' degrees of freedom. In order to extract monopoles, we divide plaquette angles 
6^ into two parts; physical fluxes (9^„ G (— it, it] and Dirac strings lim^y. 

Onu = <V + Zim^ (3) 

n^u £ [0, ±1,±2] is integer-valued and corresponds to the number of the Dirac strings pen- 
etrating the plaquette. We can now define the integer-valued DeGrand-Toussaint monopole 
current m^{x) [S] in a gauge- invariant manner; 

m^x) = -e^ VK xA+n K x(x), (4) 

where is a forward and backward derivative operator on a lattice, respectively. The 
monopole currents m^x) satisfy the conservation law A~77i M (x) = and hence form closed 
loops. 

Monopole currents can be unambiguously classified into monopole clusters C l man . In 
Ref. [2H]j it was conjectured that the largest monopole cluster C^ on occupies most of monopole 
currents and only the largest monopole cluster is relevant for the color confinement. In Fig.[TJ 
we show the scattered plots of the numbers of the sites occupied by the i-th largest monopole 
cluster at each (3, which are obtained with 48 configurations at each (3. At /3=0.99 and 1.01, 
we find prominently large monopole clusters of the length of about 10000, whereas we find 
no large cluster at 0=1.03. The second largest cluster C^ on is much smaller than the largest 
cluster Cj^ on in the confinement phases, especially at 0=0.99, which is the same tendency as 
that reported in Ref. [28]. This feature may indicate the possible interrelation in compact 
QED between the nonperturbative feature and the appearance of large monopole clusters 
which cover almost entire volume. 
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Figure 1: The numbers of the sites occupied by the i-th largest monopole cluster are plotted, 
which are obtained with 48 configurations at each (3. 



2.2 Overlap fermion 



For the fermion action, we employ the overlap formalism 
D is constructed as 



] . The overlap-Dirac operator 



D = p[l + j 5 sgn(H w )] = p 



f + 75 



H 



w 



(5) 



and realizes the exact chiral symmetry on a lattice satisfying the Ginsparg- Wilson rela- 
tion [SUES] 

l5 D + D l5 =p- 1 Dj 5 D. (6) 

Here, Hw = 75 (Dw — p) is the hermitian Wilson-Dirac operator defined with the standard 
Wilson-Dirac operator Dw- The "negative mass" p is chosen in the range of < p < 2, 
which we set f .6 throughout this paper. We approximate the sign function sgn(Hw) by f50 
degrees' Chebyshev polynomial [33], treating O(200) lowest eigenmodes of Hw exactly. In 
the present analysis, we impose the periodic boundary conditions in all the spatial direc- 
tion for the fermion fields, whereas the anti-periodic boundary condition is imposed in the 
temporal direction. We compute lowest 50 eigenpairs for each f3 implementing a restarted 
Arnoldi method [34] ■ All the eigenvalues Ai a t of D, which lie on a circle with the radius of 
p in a complex plain, are stereographically projected onto the imaginary axis via Mobius 
transformation [35], 



A 



1 - A lat /2p' 



(7) 



3 Low-lying Dirac modes 

We briefly survey the properties of Dirac eigenmodes in this section. 

We first show the histograms H CV (X) of Dirac eigenvalues in Fig. [2] where the horizon- 
tal axis denotes — iX. Since a non-zero eigenvalue is always accompanied by its complex 
conjugate, the figure is symmetric about the H ev axis. We find exact zero modes in the 
confinement phase at /3=0.99 and l.Of , but on the other hand no zero mode is seen in the 
Coulomb phase at (3—1.03, which is consistent with the previous works ^23, 25 . The numbers 
of zero modes found in 48 gauge configurations at each (3 are listed in Table. [T] The other 
prominent difference can be seen in the spectral density at the spectral origin. At the strong 
coupling (/3=0.99 and 1.01), the density p ev {\) of near-zero modes is rather dense, while it 
rapidly goes to zero at /3=1.03 as ~ |A CV | 3 - These non- vanishing eigenvalue densities at A ~ 
are directly connected to non-zero chiral condensate via Banks-Casher relation [4] as 

(^)ocpev(O). (8) 

Such dense spectral density is then considered as a signal of the broken chiral symmetry in 
the strong-coupling compact QED. 
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Figure 2: The histograms H ev (X) of the eigenvalues A of the overlap-Dirac operator D are plotted 
for each (5. The horizontal axis denotes — iX. All the eigenvalues Aj a t lying on a circle in a complex 
plain are stereographically projected onto the imaginary axis via Mobius transformation, A = 

l-A lat /2p- 
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Table 1: The numbers of exact zero-modes found in 48 gauge configurations are listed. The i-th 
column gives the number of configurations with 0, 1, 2 zero-mode(s), respectively. 



We next extract the inverse participation ratio (IPR) for each eigenmode at each (3. The 
IPR I(X) is denned as 

I(X) = vY / Pipa(x) 2 , pipr(x) =^|Va(*)| 2 . (9) 

x a, a 

Here, V denotes the system volume and ip\ (x) is the eigenfunction associated with an eigen- 
value A normalized as \ip\(x)\ 2 = 1. The Roman and Greek alphabets a and a are the 
indices for a color and a spinor, respectively. The density Pipr(cc) is obtained by locally 
summing up the absolute square of each component of an eigenfunction ip\(x) only over its 
color and spinor indices. The IPR is unity when ip\{x) maximally spreads over the system 
and equals to V in the case when "4>\{x) lives only on a single site, reflecting the spatial distri- 
bution of the eigenfunction tp\(x). As one can see in Fig.[3l zero modes exhibit much larger 
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Figure 3: The scattered plots of inverse participation ratios of low-lying Dirac modes at /3=0.99, 
1.01 and 1.03. The horizontal axis denotes —iX, the associated eigenvalue to each Dirac mode. 

IPRs than the other low- lying modes (near-zero modes) . This tendency implies the stronger 
localization property peculiar to zero modes in compact QED. The IPRs for near-zero modes 
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have /3-dependence; the IPRs at larger (3 are smaller than those at smaller (3 as a whole, 
which implies that near-zero modes get delocalized as we increase (3. This derealization 
property of near-zero modes was also found in Ref. [25) . Dirac modes in the weak coupling 
limit are expected to be plain waves and to be completely delocalized. The p(A) and the 
IPR at /3=1.03 supports this expectation. 

4 Histograms of Dirac eigenmodes 

Our main concern in this paper is a relationship between Dirac modes and monopoles. 
Though the IPR surely reflects one aspect of a spatial distribution of an eigenfunction, it can 
say almost nothing about the shape of an eigenmode. The difficulty in such an analysis comes 
from the difficulty in a quantitative evaluation of the detail of eigenfunction distributions. 

We first concentrate ourselves on the eigenfunction-density histograms H^ y (p^). The 
density p^ = pipr here is merely the same quantity as the IPR density Pipr of a Dirac 
eigenfunction ip\(x). The open bars in Fig. [4] show the typical density histograms H^ 1 
counted at all the sites for both zero and near-zero modes at ,3=0.99. The histograms H^ 1 
for zero modes have a relatively long tail at large p^p and show a large height at ~ in 
contrast to those of near-zero modes, which also exhibits the stronger localization property 
of zero modes. 

We next investigate another type of histogram H™ on (p^), which is defined by counting 
only u Pl p on monopoles". H™ on is here normalized as / °° H™ on (t)dt = f™ Hf(i)dt. Taking 
into account that monopoles are by definition located at the intermediate points x = x + 
(1/2, 1/2, 1/2, 1/2), p^'s are simply redefined by averaging them around intermediate points 

x as 

p 4 ,^ Pi ,{x)= Pip( x )/ 24 ( 10 ) 

\x— x\ — 1 
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Figure 4: The density histograms H^ 1 and H™ n at /3=0.99 are plotted as the function of p^. 
The histograms H^ 1 , which are counted over all the site, are shown as open bars. The histograms 
H™ 11 , which are counted only on monopole currents are shown as solid lines. 

In Fig. 21 we show the histograms -ff™ on for both zero and near-zero modes counted on 
the monopole world lines at (3—0.99 as solid lines. At small p^,, we can find a remarkable 
difference especially for near-zero modes. At p$ < 0.00005, the histogram H™ on counted 
only on monopoles takes a larger value than H^ 1 counted over all the sites. On the other 
hand, at p^ > 0.00005, H™ on tends to be smaller than H^ 1 . This tendency indicates that 
monopoles "run" on the sites where the density p^ of Dirac eigenfunctions is small. Or we 
can say that near-zero Dirac modes are localized avoiding monopoles. 

We define and investigate the histogram ratios R-^(p^) in order to evaluate the correla- 
tions in a semi-quantitative way; 
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This quantity R 1 p{p 1 p) equals to 1, if there is no correlation between the spatial fluctuations 
of Dirac modes and monopoles. In the case when a positive (negative) correlation exists 
between the spatial fluctuations of Dirac modes and monopoles, R^,(p^,) > 1 at smaller 
(larger) p^ and R^ip^,) < 1 at large (smaller) p^ hold. 




Figure 5: The histogram ratios R^(p^p) for zero (upper two panels) and near-zero modes (lower 
three panels) at /3=0.99, 1.01 and 1.03 are plotted in the range of < p^ < 0.0002. We draw a 
line at R^(p^) = 1 for reference. 

We here evaluate R^ only with p^, < 0.0002, because histograms are quite sparse 
at p^p > 0.0002 and this simple analysis cannot be reliable. (For the R^,(p^) for near-zero 
modes, we simply average R over all the available near-zero modes.) The ratios (to 1) 
of the squared norms of the Dirac modes obtained by summing up only p^ > 0.0002 are 
12(7)%, 15(11)% for zero modes at /3=0.99, 1.01, and 1.4(3)%, 0.5(3)%, 0.1(1)% for near- 
zero modes at /3=0.99, 1.01, 1.03, respectively. As expected from the localization property of 
zero modes, their ratios are larger than those of near-zero modes. In spite of the truncation 
at p^ ~ 0.0002, the analysis in the range of p^ < 0.0002 might be enough to grasp the nature 
of Dirac modes, at least for near-zero modes, since p^ > 0.0002 are responsible only for a 
few % of the total squared norms. Fig. [5] show the histogram ratios R^(p^p) plotted in the 
range of p^ < 0.0002, respectively, for both zero and near-zero modes at /3=0.99, 1.01 and 
1.03. 



4.1 Results for zero Dirac modes 

The exact zero Dirac modes in compact QED have been studied with regard to their relation 
to topological excitations [23] , and to their spatial distributions [25 . The exact zero modes 
in compact QED can be found only in the confinement phase [23l[25]. No detailed correlation 
between any of the topological phenomena and the zero-mode degeneracy of the overlap-Dirac 
operator was found in Ref. ,23]. 

While we have only a few dozen of zero modes in 48 gauge configurations and the statis- 
tical errors are consequently rather large, the spatial fluctuations of zero modes seem to have 
no remarkable correlation with monopolc currents. There are naively three possibilities: 

1) The eigenfunctions fluctuate in a very complicated way forming multifractal structures. 

2) The zero-mode's fluctuation (especially p^ < 0.0002) simply have little correlation with 
monopoles and would be governed by some other possible objects. 3)The apparent correla- 
tions between zero modes and monopoles appear only in the sector of p^, > 0.0002, which 
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we have discarded. 

Even if the relationship between monopoles and zero modes is originally simple, the 
complicated monopole structure could lead to the complex shape of zero modes, which our 
present simple analysis may not detect. (The monopole currents in fact occupy almost half 
of the total volume of 12 3 x 12 = 20736 and their structures themselves are complicated.) 
Moreover, the averaging in Eq. [TU]can be harmful in this case. In order to verify this possi- 
bility 1), it would be needed to compute multifractal dimensions of the modes by analyzing 
the scaling of eigenfunction moments, which is left for further study. In the case of 2) and 
3), our present analysis can cast no more light on the nature of zero modes. The authors in 
Ref. [36] conjectured that Dirac zero modes mainly appear at the intersections of vortexes. 
In such a rather complicated case, we need another analysis. 

In any case, we would need more statistics and detailed analyses in order to make definite 
conjectures upon zero modes. 

4.2 Results for near-zero Dirac modes 

As for near-zero Dirac modes, one can see the outstanding anti- correlation in Fig. [5] The 
histogram ratio R^,(p^) is clearly larger than 1 at small p^ and less than 1 at large p^ 
beyond the statistical errors. Although the values of R^(p^) themselves are slightly different 
among three /3's, /3=0.99, 1.01 and 1.03, the bulk properties are surely the same, which 
means the existence of the universal anti-correlations between near-zero Dirac modes and 
monopole currents around the critical coupling f3 c . The important point is the universality 
of the anti-correlation between near-zero modes and monopoles. This simple rule does not 
drastically change before/after the phase transition. Near-zero modes are "scattered" by 
monopole clusters. 



5 Discussions and Speculations 

5.1 Dirac eigenvalues and Chiral condensate 

As we have seen in Sec. [3l Dirac eigenvalues have non- vanishing density p cv (0) at the spectral 
origin in the chiral broken phase, and this non-vanishing value is directly connected to the 
non-zero chiral condensate (tpip) via the Banks-Casher relation. The non- vanishing p cv (0) is 
generated by the "repulsive force" among the eigenvalues: The eigenvalues of the overlap- 
Dirac operator lying on the circumference of the circle in a complex plain repel each other, 
and consequently near-zero eigenvalues are pushed towards the spectral origin forming the 
non-zero p C v(0). Namely, the driving force of the chiral symmetry breaking is the "repulsive 
force" among Dirac eigenvalues. On the other hand, when p ev (0) is zero, the repulsive force 
is considered to be weaker or zero. The magnitude of this force can be revealed by level 
statistics. 




012340123401234 
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Figure 6: The unfolded nearest-neighbor level spacing distributions fl a t(s) at j3 =0.99, 1.01 and 
1.03. The solid lines denote the Wigner distribution function Pwi s (s) = ||s 2 exp(— ^s 2 ) and 
the dashed lines the Poisson distribution function Pp Q i(s) = exp(— s). 
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We show in Fig. [S] the unfolded nearest-neighbor level spacing distributions P\ a t(s) of the 
low-lying Dirac eigenvalues obtained at (3=0.99, 1.01 and 1.03. The solid lines denote the 
Wigner distribution function Pwig(s) = p-s 2 exp(— |:S 2 ) and the dashed lines the Poisson 
distribution function Pp i(s) = exp(— s). The Wigner distribution function -Pwig(s) is a good 
approximation of the original distribution function obtained by the random matrix theory 
(RMT) for the chiral unitary ensemble (chUE). The Poisson distribution Pp i(s) appears, 
for example, in the system where eigenenergy levels have no correlation with each other. 

The level spacing distributions P\ a t(s) at ,3=0.99 and 1.01 show the good coincidence 
with Pwig(s)- The manifestation of the Wigner distribution, at (3=0.99 and 1.01 in the 
chiral broken phase, implies the level repulsion among the low-lying eigenvalues, which is 
consistent with the non- vanishing eigenvalue density p cv (0) at the spectral origin. On the 
other hand, at /3=1.03 in the chiral restored phase, Pi a t(s) coincides with neither the Wigner 
distribution Pwig(s) nor the Poisson distribution Pp os (s). P\ a t(s) at (3=1.03 seems to be on 
the way from the Wigner to the Poisson distribution. As a remarkable fact, you can find the 
long tail of P\ & t(s) stretching even at s = 4. This tendency seen at (3=1.03 shows a weaker 
repulsion among the low-lying eigenvalues than at (3=0.99 and 1.01, which leads to small or 
zero eigenvalue density p e v(0) at the spectral origin. The Wigner-Poisson transition in the 
nearest-neighbor level spacing of the low- lying Dirac eigenvalues appears at (3 > 1.01, almost 
at the same time as the confinement-deconfinement phase transition. 

So far, we have mainly shown the numerical facts obtained in the present analysis. In 
the next subsection, we make a possible speculations about the chiral transition mechanism 
in compact QED. 

5.2 Possible speculations 

What is the origin of a Wigner-Poisson transition in the neighboring level spacing distribu- 
tions of the low-lying Dirac eigenmodes, which is responsible for the chiral phase transition? 
The Wigner/Poisson distributions in neighboring level spacing distributions can be also found 
in classically chaotic/regular systems. Berry and Tabor [37] drew the conclusion that in a 
classically integrable system with more than one degree of freedom the level spacing dis- 
tribution of quantum spectra obeys the Poisson distribution. We can also find in Ref. [38] 
the famous conjecture by Bohigas, Giannoni and Schmit stating that the level statistics of a 
quantum system whose classical counterpart is chaotic do not show a dependence on the de- 
tails of the dynamics but depend only on the global symmetry of the system. Level statistics 
in such a system coincide with those obtained by the random matrix theory with the same 
global symmetry. (This conjecture is not always true. A quantum kicked rotor does not 
satisfy the conjecture.) The chaoticity or the complexity of the system leads to the Wigner 
distribution in level spacing. 

The key player for the Wigner-Poisson transition in compact QED would be monopoles. 
As was seen in Sec. [4] there exist universal anti-correlations between monopoles and near- 
zero Dirac eigenmodes. In the system where large monopole clusters exist, the near-zero 
modes are "scattered" by monopoles in a very complicated way, which results in the Wigner 
distribution of the neighboring level spacing. When there exist only small monopole clusters 
in a system, this complexity is much weakened and the level spacing exhibits the Poisson- 
like distribution. The level dynamics of low-lying Dirac eigenvalues are simply controlled by 
monopoles' clustering-declustering feature. 

Then the fact that the chiral symmetry breaking is always accompanied by the charge 
confinement in compact QED could be naturally understood from the microscopic view- 
point. On one hand, the monopoles' clustering-declustering feature is surely responsible for 
the confinement-deconfinement transition. On the other hand, the monopoles' clustering- 
declustering transition gives rise to the change in the complexity of the system and conse- 
quently in the level dynamics of low-lying Dirac eigenmodes, which leads to the change in 
the eigenvalue density at the spectral origin, or equivalently in the chiral condensate i'tptp). 

We finally note that the chaoticity and the level spacing distribution in gauge theories 
were also investigated by several groups [221 13"9l |4"01 HTl |4"21 |4"31 EH] . 
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6 Summary 



We have studied the properties of low-lying Dirac modes in quenched compact QED at 
/3=0.99, 1.01 and 1.03, employing 12 3 x 12 lattices and the overlap formalism for the fermion. 
We have found several features worth noting: 

• The nearest-neighbor level spacing distribution of the Dirac eigenvalues coincides with 
the Wigner distribution at (3 = 0.99, 1.01 < (3 C , and it seems to be on the way from 
the Wigner to the Poisson distribution at /3 = 1.03 > (3 C , which is consistent with the 
non- vanishing (vanishing) chiral condensate at (3 < (3 C (/? > f3 c ). 

• The chiral phase transition in compact QED is very different from that in QCD at finite 
temperature. The near-zero Dirac modes are much delocalized in the chiral restored 
phase in compact QED, whereas they are strongly localized in finite temperature QCD 
exhibiting the Anderson-transition(AT) of the vacuum [45"1 146] . 

• Near-zero modes have been found to have universal anti-correlations with monopole 
world lines below/above the critical (3. 

The anti-correlation between monopoles and near-zero modes indicates that near-zero 
modes are "scattered" by monopole currents in a complicated way, which is considered as 
the origin of the critical behaviors of low-lying Dirac eigenmodes. The chiral phase transition 
in 4D compact QED might be then associated with the complexity of the vacuum, which is 
brought about by the clustering-declustering transition of monopoles. 
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